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Abstract 

A new linked cluster expansion for the calculation of ground state observables of complex nuclei 



I; using the V8' 



potential the ground state 



with realistic interactions has been developed ODD 
[~^. , energy, density and momentum distribution of complex nuclei have been calculated and found to 

o 

be in good agreement with the results of obtained within the Fermi Hyper Netted Chain, and 
in ■ Variational Monte Carlo 6] approaches. Using the same cluster expansion, with wave function 

and correlations parameters fixed from the calculation of the ground-state observables, the semi- 
inclusive reaction of type A(e, e'p)X has been calculated taking final state interaction effects into 
account within a Glauber type calculation as in Ref. [_|]; the comparison between the resulting 
distorted and undistorted momentum distributions provides an estimate of the transparency of 
the nuclear medium to the propagation of the__h.it proton. The effect of color transparency has 



also been considered within the approach of |9j; it is shown that at high values of Q 2 finite 
formation time effects strongly reduce the final state interaction, consistently with the idea of a 
reduced interaction of the hadron produced inside the nucleus |lo| . The total neutron-nucleus cross 
section at high energies has also been calculated \l\ by considering the effects of nucleon-nucleon 
correlations, which are found to increase the cross section by about 10% in disagreement with the 



experimental data. The inclusion of inelastic shadowing effects of Refs. 



12 



3 



decreases back the 



cross section, leading to a good agreement between experimental data and theoretical calculations. 
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I. INTRODUCTION 



The knowledge of the nuclear wave function, in particular its most interesting and poorly 
known part, viz the correlated one, is not only a prerequisite for understanding the details 
of bound hadronic systems, but is becoming at present a necessary ingredient for a correct 
description of medium and high energy scattering processes off nuclear targets; these in fact 
represent a current way of investigating short range effects in nuclei as well as those QCD 
effects (e.g. color transparency, hadronization, dense hadronic matter, etc) which manifest 
themselves in the nuclear medium. The necessity to treat nuclear effects in medium and 
high energy scattering within a realistic many body description, becomes therefore clear. 
The problem is not trivial, for one has first to solve the many body problem and then has 
to find a way to apply it to scattering processes. The difficulty mainly arises because even 
if a reliable and manageable many-body description of the ground state is developed, the 
problem remains of the calculation of the final state. In the case of few-body systems, a 
consistent treatment of Initial State Correlations (ISC) and final state interaction (FSI) is 
nowadays possible at low energies by solving the Schrodinger equation for the bound and 
continuum states but, at high energies, the Schrodinger approach becomes impractical and 
other methods have to be employed. In the case of complex nuclei, much remains to be 
done, also in view that the results of very sophisticated calculations (e.g. the variational 
Monte Carlo ones [(|), show that the wave function which minimizes the expectation value 
of the Hamiltonian provides a very poor nuclear density; moreover, the structure of the best 
trial wave function is so complicated, that its use in the calculation of various processes 
at intermediate and high energies appears to be not easy task. It is for this reason that 
the evaluation of nuclear effects in medium and high energy scattering processes is usually 
carried out within simplified models of nuclear structure. As a matter of fact, initial state 
correlations (ISC) are often introduced by a procedure which has little to recommend itself, 
namely the expectation value of the transition operator is evaluated with shell model (SM) 
uncorrelated wave functions and the final state two-nucleon SM wave function is replaced 
by a phenomenological correlated wave function; to date, a consistent treatment of both 
ISC and FSI in intermediate and high energy scattering off complex nuclei is far from being 
completed, so that a quantitative and unambiguous evaluation of the role played ISC is 
still lacking. For such a reason we have undertaken the calculation of the ground-state 
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properties (energies, densities and momentum distributions) of complex nuclei within a 
framework which can be easily generalized to the treatment of various scattering processes, 
keeping the basic features of ISC as predicted by the structure of realistic Nucleon-Nucleon 
(NN) interactions. This paper is organized as follows; in Section [HI the cluster expansion 
method is introduced, and the relevant ground state properties of 16 and 40 C a calculated. 
In Section UTTl the semi-inclusive A(e, e'p)X reaction off complex nuclei targets is considered, 
and the FSI is calculated within the Glauber and Finite Formation Time pictures, taking 
advantage of the wave functions obtained in Section |H] In Section IIVI the total neutron- 
nucleus cross section is calculated with the same wave functions, and the result of taking 
into account correlations and inelastic shadowing effects are discussed. 



II. THE CLUSTER EXPANSION 



We write the nuclear Hamiltonian in the usual form, i.e.: 

H = f + V = E V, 2 + E v 2 (x hXj ) , (1) 

ZM N i= i i<j 

where the vector x denotes the set of nucleonic degrees of freedom, and the two-body 
potential J2i<j ^2{ x ii x j) is of the form 

N 



v = E^y<f , (2) 

n=l 



where = — r*^| is the relative distance of nucleons % and j, and n, ranging up to iV = 18, 
labels the state-dependent operator O^ 1 - 



(n) 



1 , (Ti • o-j , Sij , (5 • JD)„ , L 2 , .. <g) [1 , n ■ tj] . (3) 



The evaluation of the expectation value of the nuclear Hamiltonian (pQ) is object of intensive 
activity which, in the last few years |~L4j |. has produced considerable results; nevertheless, 
the level of complexity of the obtained wave functions is such that they cannot be used in 
scattering problems with reasonable ease. Our goal is to present a more economical, but 
effective method for the calculation of the expectation value of any quantum mechanical 
operator A 

{A) ~ (1>o\1>o) ' () 



with ip having the following structure 

ip (x 1 ,...,x A ) = F(xi, ...,x A ) <j> {xi, ■~,xa) , (5) 

where O is a Shell-Model (SM), mean-field wave function, and F is a symmetrized correlation 
operator, which generates correlations into the mean field wave function. According to the 
two-body interaction of Eq. (J2J), the correlation operator is cast in the following form: 

A 

F(xi, x 2 ... Xa) = S]Jffa) (6) 

i<j 

with 

JV 

ffa) = E f (n) fa) f {n) fa) = f W fa) 6ff . (7) 

n=l 

In the present paper we are going to introduce a cluster expansion technique in order to 
evaluate Eq. (J3J); the solution can be found by applying the variational principle, with the 
variational parameters contained both in the correlation functions and in the mean field 
single particle wave functions. The expectation value (A) defined in (J1J) can be expanded 
in the framework of the cluster expansion developed in Refs. Q, Q Q and, at first order, it 
reads as follows 

(^)l = (to I E (fijAfii - A ) I <t>o) - (A) o (<t> I E (4-4- - l) I 0o> , (8) 

where (A) Q is given by (0 O \A\ </>„) The 2nc? order term can straightforwardly be obtained by 
the same technique used to derive Eq. (jHJ) Given the two-body interaction as in Eq. (J2J), the 
expectation value of the Hamiltonian can be written in the following way: 

h 2 



ri=r{ + 



+ Y, f dr 1 dr 2 v^(r 12 )p^ ) (r 1 ,r 2 ), (9) 

n J 

where p^(ri,r[) and p^(r 1 ,r 2 ) are the non-diagonal One Body Density Matrix (OBDM) 
and the (spin and isospin dependent-; see Ref. 1[) Two Body Density (TBD) matrices, 
respectively These can be calculated from the ground state wave function by inserting in 
Eq. (@J) the corresponding operators. The knowledge of the OBDM and TBD matrices 
allows one to calculate, besides the ground-state energy, other relevant quantities like e.g. 
the density distribution: 

p(r) = pU(r l = r' 1 = r), (10) 



and the nucleoli momentum distribution, defined as: 

n(k) 



(27T) 



^3 Jdr^e^-^p^r',). 



(11) 



The results of calculation of the ground state energy, the charge density and the two-body 
density and momentum distribution using the realistic V8' interaction |4j is discussed in 
detail in Ref. Q, 0, jj. In Fig. El we show the effect of correlations on the charge density 
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FIG. 1: The charge densities of ie O (left) and 40 Ca (right) calculated within the cluster expan- 
sion compared with the experimental data. The mean field result and the results of the FHNC 
calculation of Ref. [|| are also shown. 

and the two-body density of 16 0. In the case of charge density, we have shown the three 
contributions which arises from the first-order cluster expansion, i.e., the mean field, hole 
and spectator contributions, which are shown in Fig. 0] in a diagrammatic picture. In the 
case of two-body density, we compare the mean field result with the ones obtained taking 
into account the only central correlation (f\ approximation) and the first six correlations 
(N = 6 in Eq. (jZJ); / 6 approximation) . 



III. APPLICATIONS. I - THE SEMI-INLCUSIVE PROCESS A(e,e'p)X 

Using the cluster expansion developed in in the Section InJ we have calculated the semi- 
inclusive A(e,e'p)X process in which an electron with 4-momentum k\ = {ki,iei}, is scat- 
tered off a nucleus with 4-momentum Pa = {0, iMa} to a state ki = {k.2, ^2} and is detected 
in coincidence with a proton p with 4-momentum p = {p, iE p }; the final (^4 — 1) nuclear sys- 
tem with 4-momentum P x = {Px,iE x } is undetected. The cross section for the exclusive 
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process A(e, e'p)B can be written as follows 

da 



Ka ep P D (E m ,p m ) (12) 



dQ 2 dvdp 

where K is a kinematical factor, a ep the off-shell electron-nucleon cross section, and Q 2 = 
|q| 2 — v 2 the four momentum transfer. The quantity P£,(E m ,p m ) is the distorted nucleon 
spectral function which depends upon the observable missing momentum p m = q — p 
(Pm = ~~ k when the FSI is absent) and missing energy E m = v — T v — Ta-i- In the semi- 
inclusive A(e, e'p)X process, the cross section (j!2)) is integrated over the missing energy 
E m , at fixed value of p m and becomes directly proportional to the distorted momentum 
distribution 

n D ( Pm ) = (2tt)- 3 j e i ^'^-^p D (r 1 ,r' 1 )dr 1 dr' 1 (13) 

where 

Mri,ri) - W^) (14) 

is the distorted one-body mixed density matrix, S is the S-matrix describing FSI and 

Pi(ri,r[) = E *(r, - f *(rj - fi) II - ^i) ( 15 ) 




FIG. 2: The momentum distributions of 16 (left) and 40 Ca (right) calculated within the cluster 
expansion, compared with the mean field result and with the FHNC approach of Ref. -5] and the 
VMC approach of Ref. j|. 
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FIG. 3: Effect of correlations versus the non-correlated result for the charge density (left) and the 
two-body density (right) of ie O calculated within the cluster expansion. 

is the one-body density matrix operator; the primed quantities have to be evaluated at 
with i = 1, A. The integral of nc(p m ) gives the integrated nuclear transparency T 




FIG. 4: Diagrammatic representation of the one body mixed density matrix p(ri,r'i) in the lowest 
order of the //-expansion (Eq. (J10j) b a) is the shell model contribution and b) and c) the hole and 
spectator contributions, respectively. The direct and exchange contributions are shown on the left 
and right sides of the Figure, respectively. 
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where Pd{?) = Pc( r i — r 'i = r ) and AT originates from the FSI. In Ref. |7[ Eq. (fTBj) has 
been evaluated using a Glauber representation for the scattering matrix S, viz 

A A 

S - S G ( ri . . . r A ) = I] C(r 1; r,-) = nil - % - 2i)r(bx - b,-)] (17) 

where bj and 2j are the transverse and the longitudinal components of the nucleon coordinate 
Yj = (hj,Zj), r(h) the Glauber profile function for elastic proton nucleon scattering, and 
the function 6(zj — z\) takes care of the fact that the struck proton "1" propagates along a 
straight-path trajectory so that it interacts with nucleon "j" only if Zj > Z\. Generalizing 
the same cluster expansion described in Section [H] to take into account Glauber rescattering, 
we have obtained the distorted nucleon momentum distributions nz>(p m ) = noipm, where 
9 is the angle between q and p m ; the results for 16 and A0 Ca are presented in Fig. 




FIG. 5: The distorted momentum distribution, n.£>(p m ) = nD(p m ,6) (8 = qp m ), of ie O and 40 Ca, 
obtained by Eq. 

The Glauber multiple scattering picture can be implemented by taking into account Finite 
Formation Time (FFT) introduced in Ref. 0, where it has been shown that at the values 
of the Bjorken scaling variable x = Q 2 jlmv ~ 1, FFT effects can be treated in a simple 
way, i.e. by replacing the Glauber operator (Eq. (JT7j) ) with 

A 

S FFT (r u t a ) = n (l - J(*i ~ %)r(b a - b,)) , (18) 

3=2 

where 

1 - e z ^^) (19) 



m being the nucleoli mass and M 2 = m* 2 — m 2 is a parameter describing the average 
excitation energy of the ejectile. It can be seen that at sufficiently high values of Q 2 , J — ► 
and the FSI vanishes; this new mechanism for the description of FSI is inspired by the QCD 
prediction of color transparency, which states that a pointlike, color-less particle, such as 
the object singled out by the virtual photon interaction with the nuclear medium at high 
Q 2 , should have reduced cross-section with the nuclear medium as long as it does not evolve 
into a physical proton inside the nucleus. We have calculated the effect of FFT on the 
forward-backward asymmetry, defined in Ref. as 

A n D {p m ,e = 0°) - n D ( Prn ,9 = 180°) 

FB n D (p m , 9 = 0°) + n D {p m , 9 = 180°) 1 ' 

and the results are shown in Fig. El 

IV. APPLICATIONS. II - THE TOTAL NEUTRON-NUCLEUS CROSS SECTION 



The neutron-Nucleus (n — A) total cross section, a l °X is defined by the optical theorem as 

4-7T 

a tot = T Im[F 00 (0)} ; (21) 
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FIG. 6: The forward-backward asymmetry (Eq. I2()|) for the process A(e, e'p)X, with A = 16 and 
A = 40, calculated within the Glauber approach (solid thick line) taking into account the Finite 
Formation Time (FFT) mechanism of various values of the 4-momentum transfer. 
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where, within the Glauber eikonal approximation, the elastic scattering amplitude has the 
following form: 



ik 



F 00 (q) = — I e^^ \{l-S G }\^ ) 



ik 
2^ 



ci 2 & n e iqbn [l - e ix °p t(bn) 



with 



^ iXopt(bn) = / n^i^-i) Mri,---,TA)\ 2 s(j y £ l T j 



(22) 



(23) 



where Vj = (sj, Zj) and b n is the neutron impact parameter. As it is well known the squared 
wave function can be written in terms of density distributions as follows jll|: 



A A A 

^ ( ri ,...,r A )f = J] Pfa) + £ A(ri,rj) J[ p{r k ) + 

3=1 «<i fe^(y) 



(24) 



where p(r) is the one-body density distribution and the two-body contraction A is defined 
in terms of one- and two- body density distributions, p^ir^r^): 



A(n,r a ) = p (2) (r 1; r 2 ) - pW(n)pW(r 



(25) 



Usually Glauber-type calculations are based upon the single density approximation, consist- 
ing in disregarding all terms of Eq. (}2*4"|) but the first one. We have considered, for the first 
time, also the second term of the expansion ()24j) i.e. the effects of two-body correlations. In 



llj ; whereas for 



the case of 4 if e we have also considered three- and for-body correlations 
4 He Eq. (}24j) has been used, for heavier nuclei we have used its optical limit (A » 1) in 
Eq. (|2SI) i.e. 



Ax op t(b n ) 



exp 
+ A 



A J dr lP ( ri )T(b n 



si) + 



2 / dr x dr 2 A(r 1; r 2 ) T(b n - s x ) T(b n - s 2 ) 



(26) 



l-Jdr lP (r 1 )T(b n -s l ) 

which already for A > 12 reproduces the results based upon Eq. almost exactly. We 
have evaluated the one- and two-body density matrices appearing in Eq. (pjj) starting from 
the realistic wave functions obtained in Section ITT1 

The results of calculations are shown in Fig. |S1 It can be seen that the inclusion of 
correlations in the target wave function produce an enhancement of the cross section of 
about 10% with respect to the mean field result, increasing the disagreement with the 
experimental data. It is well known however that at high energies diffractive scattering of 
the projectile, depicted if Fig. [7[ plays a relevant role. We have evaluated such an effect 
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a) 



FIG. 7: The neutron- nucleus total cross section diagrams, a) Glauber double scattering; b) inelastic 
shadowing effects. 



(also known as Gribov inelastic shadowing) according to Refs. |l2t Il3{]. It can be seen that 
inelastic shadowing effects play indeed a relevant role to bring theoretical calculations in 
good agreement with experimental data. 
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FIG. 8: The double scattering diagrams in the neutron-nucleus total cross section. 
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V. SUMMARY 



We have developed a method which can be used to calculate scattering processes at 
medium and high energies within a realistic and parameter-free description of nuclear struc- 
ture. Our calculations followed the following strategy: i) the values of the parameters 
pertaining to the correlation functions and the mean field wave functions, have been ob- 
tained from the calculation of the ground-state energy, radius and density of the nucleus, to 
a given order of the expansion, using realistic nucleon-nucleon interactions; ii) using these 
parameters we have calculated, by a proper generalization of the cluster expansion, the 
distorted momentum distributions, the nuclear transparency and the total neutron-nucleus 
cross section. The method we have developed appears to be a very effective, transparent 
and parameter-free one. To sum up, we have shown that, using realistic models of the 
nucleon-nucleon interaction, a proper approach based on cluster expansion techniques can 
produce reliable approximations for those diagonal and non diagonal density matrices which 
appear in various medium and high energy scattering processes off nuclei, so that the role 
of nuclear effects in these processes can be reliably estimated without using free parameters 
to be fitted to the data. 
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